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Abstract 

Randomly correlated ensembles of two quantum systems are inves- 
tigated, including average entanglement entropies and probability dis- 
tributions of Schmidt-decomposition coefficients. Maximal correlation 
is guaranteed in the limit as one system becomes infinite-dimensional. 
The reduced density operator distributions are compared with distri- 
butions induced via the Bures and Hilbert- Schmidt metrics. 
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1 Introduction 

This Letter is primarily motivated by the following question: what statis- 
tical ensemble corresponds to minimal prior knowledge about a quantum 
system? Such an ensemble may be identified as the most random ensemble 
of possible states of the system. It would provide, for example, a natural 
benchmark for assessing how "random" a given evolution process is [El]; a 
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worst-case scenario for general schemes for extracting information about the 
system [@J; and a natural unbiased measure over the set of possible states of 
the system (which would allow one to calculate, e.g., the average effective- 
ness of a general scheme for distinguishing between quantum states |J). 

For the case where the system is known to in fact be in a pure state, 
there is an obvious answer to the question. In particular, identifying minimal 
knowledge with maximal symmetry, it is natural to require that the ensemble 
be invariant under the full group of unitary transformations (thus there is no 
preferred measurement basis for extracting information). This requirement 
yields a unique probability distribution over the set of pure states of the 
system which has found applications in quantum inference Q, quantum 

chaos Q, and quantum information ^ []. 

However, as pointed out by Wootters H], there does not appear to be 
a natural generalisation of the above ensemble when the restriction of pure 
states is removed. Indeed, if general states described by density operators 
are allowed, the requirement of unitary invariance only implies that the 
probability measure over the set of possible states is a function of the density 
operator eigenvalue spectrum alone. Hence a unique probability measure can 
be specified only via some further principle or restriction, to be motivated 
on physical or conceptual grounds. 

In this Letter two possible approaches to the question are examined. 
The first is motivated by recent work of Braunstein || , and is considered in 
sections 2 and 3 below. It corresponds to assuming that the quantum system 
is randomly correlated with a second system, where the composite system is 
in a pure state. The reduced ensemble of the system is characterised by the 
distribution of Schmidt-decomposition coefficients of the composite system, 
and is explicitly calculated for the 2-dimensional case. This further allows 
calculation of the average "entanglement entropy" || of the systems. In the 
limit as the dimension of the auxilary system becomes infinite, the systems 
become maximally correlated with probability unity. 

The second approach, studied in section 4 below, is more formal in na- 
ture. It relies on choosing a metric on the space of density operators of 
the system; the random ensemble then corresponds to the (normalised) vol- 
ume element on this metric space. There are strong information-theoretic 

1 If the average density operator of the ensemble is also known, then the above- 
mentioned pure-state ensemble may be modified to give a corresponding maximally- 
random or "Scrooge" ensemble, with the property of being that ensemble of pure states on 
which measurement yields the least possible information consistent with the prior knowl- 
edge §. 
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grounds for motivating the choice of the Bures "distinguishability" metric 
jj]]. Moreover, for a 2-dimensional system, this metric corresponds to the 
conceptually satisfying case of a maximally symmetric space, with no pre- 
ferred locations or directions in the space of density operators. The ensemble 
induced by the Hilbert-Schmidt metric is also considered. 

Finally, in Section 5 comparisons are made between the above two ap- 
proaches for the two-dimensional case. It is argued that it is the second 
approach, based on the Bures metric, which yields the desired "minimal 
knowledge" ensemble in this case. 

2 Randomly correlated ensembles 

Let quantum system S, with Hilbert space Hs, be correlated with an aux- 
ilary system A, with Hilbert space Ha (for example, a similar system, a 
measuring apparatus, or the environment). A general pure state of the 
composite system then has the form 

M N 
i=\j=x 

where {| Ui)} and {| Vj)} denote orthonormal bases for Hs and Ha respec- 
tively, and M and N are the respective dimensions of Hs and Ha- 

It is always possible to choose orthonormal bases {| u*)}, {\ v*)} for Hs 
and Ha in which | ip) has the Schmidt-decomposition form |§ 

K 

\4>) = J2^\u* k )®\v* k ), (2) 

k=l 

where K <mm(M , N) , and the Schmidt coefficients {X k } are non-zero and 
unique up to permutations. These coefficients are just the (non-zero) eigen- 
values of the reduced density operators ps = tr^\ ^){ip |] and pa = 
trs[ip)(ip |] of S and A respectively, as may be verified directly from Eq.(§). 
The quantity 

E^ = - ^\ k log 2 A fc (3) 

k 

is called the "entanglement entropy" of the two systems B , and is a useful 
measure of the degree of correlation between S and A [||, y, |j| . 

It will now be assumed that S and A are randomly correlated, by which 
it is meant that the composite system is a member of the maximally random 
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pure-state ensemble discussed in the Introduction, described by a uniform 
distribution over the pure states of H$ <S> Ha- Such correlations may arise 
if the composite system is "chaotic" in the sense of Schack and Caves JLOU 
(i.e., if its state is randomised over the Hilbert space by stochastic fluctu- 
ations). They are also relevant if two observers observe the respective S 
and A components of an ensemble of pure composite systems (possibly for 
cryptographic key generation ||), in the case of minimal knowledge about 
the ensemble. Further, for M = 2, Braunstein has used such randomly 
correlated ensembles to numerically generate reduced density operators of 
S, to test the average effectiveness of a general measurement scheme for 
distinguishing between two states of S |J . 

Now, for a quantum system of dimension D, the maximally random 
pure-state ensemble over states {| a)} of the system is described by the 
probability measure |2[ Q 

D 

dn a = K D 5((a | <r> - 1) J] dRe{a d }dIm{a d }, (4) 

d=i 

where the {cr^} are the coefficients of | a) with respect to some (arbitrary) 
orthonormal basis, and the normalisation factor Kd is given by 

K D = (D- l)!/vr D . (5) 

Hence the randomly correlated ensemble is described by the corresponding 
probability measure dtl^ over the pure states | ifi) of H $ <8> Ha- 

dn^ = 2- mn K M n5({iP I V> - 1) II dPvd&j' ( 6 ) 

where the coefficients Cjj in Eq. (|l|) have the polar form (pjj) 1//2 exp(i^jj). 
It proves useful to rewrite Eq. (Q) via the definitions 

Xi = ^Pij, (7) 

j 

Pij = Pij/xi, (8) 
I cti) 



= \J Pij e*v{i(t>ij) I Vj), (9) 
j 

dn ai = 2- N K N 5({ai | ai) - VUdPijdfaj, (10) 
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where d$l ai is the uniform measure over pure states {| Oj)} of Ha- In 
particular, if for each i the variables pij in Eq. (||) are transformed to the 
variables Pn,. . .,Pi t N-i and X{, one obtains 

N 2V-1 

n = { Xi ) N - i dxi n dp^, (ii) 
j=i j=i 

where the Jacobian factor (xi) N ~ 1 is most easily evaluated by adding the 
first N — 1 rows of the Jacobian determinant to the last row. Substituting 
Eqs. (|7|)-([lTl) into Eq. (||), and multiplying by dummy terms of the form 
5((cti | a.i)— 1) dPiN, yields the final symmetric expression 

M 

dQ^ = K M n/(K n ) m 5(J2x 1 - 1) JKxtf^dxidn^ (12) 

i i=l 
for the randomly correlated ensemble. 

3 Statistical properties 

Since the measure dQ.^ is invariant under unitary transformations [Q. 0], 
the reduced ensemble of density operators ps=t r A\\ V')(V' |]> corresponding 
to system S, is similarly invariant under such transformations. It follows 
that the distribution of density operators ps is basis-independent, and hence 
that the reduced ensemble is characterised by a probability distribution over 
the eigenvalue spectrum of ps- As noted following Eq. (|2|), this spectrum is 
determined by the Schmidt-decomposition coefficients of | tp), and hence the 
corresponding probability distribution will be denoted by pm,n(M, M, ■ ■ ■)• 
Note that the symmetry between systems S and A implies that 

PM,N = PN,M- ( 13 ) 

To calculate the distribution pm,n, and hence such quantities such as 
the average entanglement entropy, note from Eqs. (jl]) and that ps 

has the general form of an M x M matrix, with coefficients 

(ui | ps | Uj) = ^/XiXj(ai | ay) (14) 

with respect to the {| U{)} basis. Hence, if a general expression for the 
eigenvalue spectrum of this matrix can be given, Pm,n(M, A2, • • •) can be 
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calculated from dCl^ in Eq. (12). This approach is successfully followed 
below for M = 2, while a less direct approach allows calculation in the limit 
as N ^ oo. 

Suppose then that M = 2. The eigenvalues of ps follow from Eq. (|i~4] ) 

as 

Ai = i(l±r),A 2 = i(l T r), (15) 

where 

r = [1 - 4xix 2 (l- | (ai | a 2 ) | 2 )] 1/2 . (16) 
As shown in the Appendix, one finds 

n (\ \ \ — (2jV-l)^(A 1+ A 2 -l) 2 ^ 

P2,at(Ai, A 2 ) - — — - _ _ — (Ai - A 2 j (AiA 2 ) , (17) 

describing a two-dimensional system S randomly correlated with an N- 
dimensional system A. 

The average entanglement entropy of S and A can be calculated from 
Eqs. (|3|) and ( |i~7| ) using standard integrals, with the final result 

_ log 2 e (2AT-1)! N ~ 2 (N-2\ (-1) s ^ 1 
\w 4iv-i (iV-2)!(JV- 1)! ^ V s / (s + 2)(2s + 3)^2t + l" 

(18) 

For the case of two randomly correlated qubits (N = 2), this yields a value of 
(log 2 e)/3~0.481 bits, which is about half of the maximum possible value of 1 
bit. In the limit as N — > oo the average entanglement entropy monotonically 
approaches this maximum (e.g., for N = 100 the average entanglement en- 
tropy is 0.99 bits). Thus maximal correlation between S and A is guaranteed 
in this limit. 

The latter result holds more generally. In fact, for arbitrary M it can 
be shown that the reduced ensemble contains only one density operator, 
M _1 1, in the limit N — > oo. The corresponding eigenvalues (and hence the 
Schmidt-decomposition coefficients) each equal M _1 , and hence from Eq. 
(||) the average entanglement entropy attains its maximum value of log 2 M, 
i.e., the systems are maximally-correlated. 

To show ps — ► M _1 l, note that integrating over the vectors {| ai)} in 
Eq. (|l~2|) yields the marginal probability distribution 

p(xi, ...,x M ) = K M n{K n )' m 5{ Xi + ... + xm- l)(xi • • • xm)*' 1 (19) 

for the diagonal elements of ps in the {| it,)} basis. Using Stirling's approxi- 
mation for n! in Eq. (0), it follows that this distribution vanishes everywhere 
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in the limit N — > oo, except for the case xt=l/M for all i. Since the reduced 
ensemble is invariant under unitary transformations, the diagonal elements 
of ps in this limit are therefore equal to 1/M relative to any basis. Choosing 
a basis in which p$ is diagonal gives ps =M _1 1 as claimed. 

The above results imply that the limit N — ► oo does not yield a par- 
ticularly "random" reduced ensemble of density operators for the system - 
indeed, it gives an ensemble with only one member. Thus, for example, the 
numerical evaluation of averages over the reduced ensemble in Section 7 of 
U for M = 2, to test the average effectiveness of a particular measurement 
scheme, is of most value for the maximally random case N = 2. 

The fact that "randomness" is in fact decreased as the dimension of 
the auxilary system increases suggests that a potential candidate for the 
minimal-knowledge ensemble discussed in the Introduction is the reduced 
ensemble corresponding to N=M, with corresponding eigenvalue distribu- 
tion pm,m (choosing iV less than M would unduly restrict the ensemble to 
density operators with M — N zero eigenvalues). However, other potential 
candidates may be generated by a second approach, as shown in the next 
Section. 



4 Metric-induced ensembles 



An ensemble of general states of a quantum system is in general de- 
scribed by a probability measure over the density operators of the system. 
Given that probability measures transform in the same way as volume ele- 
ments under co-ordinate transformations, and that volume elements are in 
general properties of metric spaces, this suggests that the distribution of 
density operators corresponding to a "minimal-knowledge" ensemble may 
be obtained from the normalised volume element induced by some natural 
metric on the space of density operators. 

A metric of particular interest is the Bures metric @ , where the infinites- 
imal distance element between two states p and p+Sp is given by 

(ds B ) 2 = 2j2(^j + Afc)" 1 I (j \5p\k) | 2 , (20) 

where p is diagonal in the orthonormal basis {| j)} with eigenvalues {Aj}. 
This metric provides a unitarily-invariant measure for distinguishing be- 
tween two quantum states, and has been strongly motivated as physically 



relevant both on measurement 12] and statistical M, 13] grounds 
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To calculate the volume element corresponding to the Bures metric, it 
is useful to decompose p+5p as an infinitesimal shift in the eigenvalues of p 
followed by an infinitesimal unitary transformation: 

p + 5p = (l + 5U)( P + 5A)(l + 5Uy 

= p + SA+[5U,p], (21) 

where (j \ 5 A \ k) = dj^dXj, and (5U)^ = —5U follows from unitarity. Note 
moreover that the infinitesimal generator 5U can generally be decomposed 
as 

SU = J2i( dx jk + idVjk) I | -h.c] (22) 

where dxjk and di/jk are real, and h.c. denotes the Hermitian conjugate of 
the expression preceding it. 

Substitution of Eqs. fnl) and <M) into Eq. (pO) yields 



(ds B f = E ^ + 4 E /^W ) 2 + (dy ]k )\ (23) 

from which one immediately extracts the volume element 

d\\...d\M -r-r . (Aj - A fc ) 2 . . 

" Kb ~ (Ai . . . a m) v2 n ^n^r*^- < 24 > 

Normalising dVs yields the desired probability distribution over the space 
of density operators. 

Since the metric is invariant under unitary transformations, the corre- 
sponding ensemble is characterised by the marginal probability distribution 
Pb(X\,. . . ,Xm) over the eigenvalue spectrum of the density operators describ- 
ing the system (see also Section 3). This distribution can be obtained from 
Eq. ( pip by integrating over the (compact) space of unitary transformations 
(parametrised by {xjk,Ujk}), and normalising, to give 

(\ x ,_ r S(X 1 + ... + X M -l) yr (A,-A fc ) 2 

Pb{Xi,. . . ,X M ) - C M 7T : — 7775 [[ > . , — , (25) 

(Ai . . . X M ) L/Z f^ k Xj + X k 

where Cm is a normalisation constant, and the condition tr[p] = 1 has been 
made explicit. 

Eq. ( ^5|) will be compared with Eq. (|l7| ) in the following Section for 
the case M = 2. This Section is concluded by noting that in principle there 
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are many possible choices of metric, each leading to a possible "random" 
ensemble. One particularly simple choice is the Hilbert-Schmidt metric, 
with infinitesimal distance element 

(ds HS ) 2 = tr[(6p) 2 }. (26) 

Following essentially the same procedure as above for the Bures metric 
(where the trace is evaluated in the {| j)} basis), one finds the corresponding 
probability distribution 

Phs(Xi, • • • , A M ) = C' M 6{\i + . . . + A M - 1) H{Xj - X k ) 2 (27) 

j<k 

for the density operator eigenvalue spectrum. This is also considered in the 
following Section for the case M =2. 



5 Two-dimensional comparisons 

The states p of a two-dimensional system may be parametrised in the Bloch 
representation as 

p=l(l + a-r), (28) 

where a is the 3-vector of Pauli matrices and r is a 3-vector of modulus 
r < 1. The eigenvalues of p are related to r as per Eq. (|i~5|). 

A distribution over p may therefore be written as a distribution over r. 
Moreover, since unitary transformations of p correspond to rotations of r, 
it follows that distributions corresponding to unitarily-invariant ensembles 
depend only on the modulus r, being uniform with respect to direction. 
Hence the distributions over r corresponding to eigenvalue distributions Eqs. 
( |j~7|) (with N = 2), ( |25|) and ( p7| ) are given respectively by 

P2, 2 (r) = 3/(4vr), (29) 
p B (r) = (AMil-r 2 )" 1 / 2 , (30) 
PHs(r) = 3/(4vr). (31) 

It is seen that the first and third distributions are uniform over the unit 
3-ball, while the distribution corresponding to the Bures metric is sharply 
peaked at the surface of the ball (corresponding to pure states of the system) . 
This raises the question of which is the more "random" ? I shall argue here 
for the latter, due to its greater symmetry. 
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In particular, the Bures metric for a two-dimensional system corresponds 
to the surface of a unit 4-ball [O], i.e., to the maximally symmetric 3- 
dimensional space of positive curvature (and may be recognised as the 
spatial part of the Robertson- Walker metric in general relativity fli|| ). This 
space is homogenous and isotropic, and hence the Bures metric does not 
distinguish a preferred location or direction in the space of density operators. 
Indeed, as well as rotational symmetry in Bloch co-ordinates (corresponding 
to unitary invariance), the metric has a further set of symmetries generated 
by the infinitesimal transformations fl4[| 

r ^r + e(l-r 2 ) 1/2 a (32) 

(where a is an arbitrary 3- vector). 

Taking the viewpoint that maximal randomness corresponds to an en- 
semble with maximal symmetry, it follows that the distribution of Eq. fl30] ) , 
in corresponding to a maximally symmetric space, is in fact more "random" 
than the distributions of Eqs. (p9|) and (31). This strongly suggests, at least 



for two-dimensional quantum systems, that the minimal-knowledge ensem- 
ble discussed in the Introduction is the one induced by the Bures metric. 

Finally, note that the existence of various candidates for the minimal- 
knowledge ensemble discussed in the Introduction begs the question as to 
whether there exists some natural physical process for generating ensembles 
of quantum systems, which can be identified with maximal randomness. 
This would allow experimental determination of the minimal-knowledge en- 
semble. This is, however, beyond the scope of this Letter. 
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Appendix 

To derive Eq. (|l7|) , note from Eqs. ([l5|) and (|l6| ) that the eigenvalue 
distribution can be calculated if the joint distribution of the variables 

X = xx,Y=\ (ai | a 2 ) | 2 (33) 



is known. From Eq. ( |l2|) the statistics of X and Y are independent, and 
hence this joint distribution has the factored form 

p(X,Y)=p(X)q(Y). (34) 



From Eqs. (12) and fl33|) one immediately has 

p(X) = K 2N (K N )- 2 [X{1 - X)} N -\ (35) 
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To find q(Y), fix an orthonormal basis {| Vj)} in Ha, and for a given unit 
vector | qi) let U be a unitary transformation which maps | a\) to | v\) and 
define | f3) =U | 02). Thus 

Y =| (ai I C/ f t/ I a 2 ) | 2 =| («i I 13) | 2 . (36) 

Since dVt a2 is invariant under unitary transformations, writing (vj \ 
P)=( Wj ) 1 / 2 exp^) yields 

d£l a2 = (Kip 

= 2~ N K N 5{(p I 0) - 1) Y[ dwjdOj 

3 

N N 

= 2~ N K N 8(J2 Wj - (1 - Y))dYdB x \{ dwjdOj 
3=2 3=2 

N N 

= 2- N K N S(J2 w j - - Y) N - 2 dYd6i J] dWjtWj 
3=2 i=2 

= 2~ 1 (^/^_i)(i - y) JV - 2 dyd0idn 7 , (37) 

where Wj =Wj/(l — Y), and <if2 7 is the uniform measure over pure states of 
the (N — l)-dimensional space spanned by {| -u 2 ). . .| vn)}- Multiplying this 
expression by d£l ai and integrating over all variables except Y then gives 

q(Y) = 7r(K N /K N ^)(l - Y) N - 2 . (38) 

Finally, substituting Eqs. (|), © and @) into Eq. (|J), the marginal 
distribution of r in Eq. ([!(]) can be calculated as 

P(r) ~ 1 r 2 fl^r" 2 (39 ) 

Plrj " 2(iV-l)!(iV-2)! r ^ 4 ) ' W 

which with Eq. (|l5| ) immediately yields Eq. (|l7[). 
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ADDITIONAL NOTES 



Equation (17) is a special case of Theorem 3, Sec. Ill of S. Lloyd and H. 
Pagels, in Ann. Phys. (NY) 188 (1988) 186. Taking M < N without loss 
of generality (see Equation (13)), one has the general formula 

PMjN (X 1 ,...,X M ) = C MjN S(J2^m-l)[l[(X m -K) 2 ][l[i X k) N ~ M ]: (40) 

m<n k 

where Cm,n is a normalisation constant. 

Equation (18) can be simplified and generalised to calculate the average 
entanglement entropy for all values of M and N, using a formula conjectured 
by Don Page and elegantly proved by S. Sen in Phys. Rev. Lett. 77 (1996) 
1-3. Again taking M < N, one has 

"! 1 M-l 

m =N+l 



As noted in the published version of the present paper (Phys. Lett. 
A 242 (1998) 123-129), a recent related preprint seeks to determine 
the "maximally noninformative" ensemble for a two-dimensional quantum 
system. This is essentially a different concept from the "maximally random" 
ensemble sought here (and identified as corresponding to the Bures volume 
measure in the 2-D case). Indeed, intuitively one would expect an ensemble 
consisting solely of the density operator proportional to the density operator 
to be least "informative", as (i) no Shannon information can be gained by 
measurement on such an ensemble; and (ii) no preferred basis can be singled 
out by measurement. 
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